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RANDOM PARKING, EUCLIDEAN FUNCTIONALS, AND RUBBER ELASTICITY 

ANTOINE GLORIA AND MATHEW D. PENROSE 



ABSTRACT. We study subadditive functions of the random parking model previously analyzed 
by the second author. In particular, we consider local functions S of subsets of R d and of 



point sets that are (almost) subadditive in their first variable. Denoting by £ the random parking 
pvj | measure in R d , and by £ R the random parking measure in the cube Qr = (— R, R) d , we show, 



under some natural assumptions on S, that there exists a constant S 6 R such that 



5 ■ lim gWglfl = lim gWgli!) = S 

fl^ + oo \Qr\ R^ + oo \Q R \ 

\j~. • almost surely. If f i-> S(Qr,() is tne counting measure of ( in Qr, then we retrieve the 

result by the second author on the existence of the jamming limit. The present work generalizes 
this result to a wide class of (almost) subadditive functions. In particular, classical Euclidean 
("V | optimization problems as well as the discrete model for rubber previously studied by Alicandro, 

CL ' Cicalese, and the first author enter this class of functions. In the case of rubber elasticity, this 

yields an approximation result for the continuous energy density associated with the discrete 
model at the thermodynamic limit, as well as a generalization to stochastic networks generated 



■i— > | 

on bounded sets. 
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1. INTRODUCTION AND INFORMAL STATEMENT OF THE RESULTS 

Renyi's model of random parking (also known as random sequential adsorption or random 
sequential packing), is defined in d dimensions as follows. A parameter po > is specified, and 
open balls -Bi,.r, B2,r, ... of radius po, arrive sequentially and uniformly at random in the d- 
dimensional cube Qr := (-R, R) d , subject to non-overlap until saturation occurs. It is known 
[PenOl] that the random parking measure £ in Qr (i.e. the point measure representing the 
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2 ANTOINE GLORIA AND MATHEW D. PENROSE 

locations of the balls at saturation) converges weakly in the sense of measures to a measure £ in 
M rf , called the random parking measure in R , and that there exists A € M + such that 

lim = A 

almost surely — which yields the existence of a deterministic averaged density of packed balls 
in M d . Perhaps of greater interest, however, is the existence of the limit of . Af , which would 
be the thermodynamic limit of the averaged density of packed balls in the domains Q R . Using 
quantitative properties of "stabilization" of the random parking measures, it is proved in [PenOl] 
that in probability, 

(1) lmi = A. 

R-^co \Q R \ 

That is, not only does the measure ( R converge weakly to £, but also the averaged density of 
packed balls converges to A (the jamming limit). Moreover, if the arrivals processes for different 
R are coupled by being all derived from a single Poisson process in space-time, the convergence 

(1) holds almost surely, as does the convergence of £ R to ( (in fact, almost surely for all r > 
we have £ R n Q r = £ n Q r for all large enough R). 

The first part of this paper is concerned with the extension of (1) to more general functions 
of £ R , besides the total measure of £ . We consider local subadditive functions S of bounded 
open sets and point sets, that is such that for all open bounded disjoint subsets D, D\, . . . , D„ 
of M d , and for every point set ( in some fixed class (i. e. satisfying the non-empty space and 
hard-core conditions, see Section 2.2), we have (see Theorem 2.3 for milder conditions) 

S(D,() = S(D,(\ D ), 

S(u? =1 A,C) < $>(A,0- 

Under the following further assumptions on S: 

• uniform boundedness: there is a constant C such that \S(D, ()\ < C\D\ for all D, (; 

• insensitivity to boundary effects: there exists < a < 1 such that iq | ~ fl = 

o(i); 

we shall prove that there exists deterministic S £ M + such that almost surely 

„, ,. S(Q R ,£ R ) S(Q R ,Q -= 

(2) lim — —— — = lim — — — : — = S. 

R-^+oo \Q R \ R-^+oo \Q R \ 

Note that (1) is a particular case of (2) for the additive function S(D, £) := ((D) = J D d((x). 
Our proof of (2) differs from the proof of (1) in [PenOl] in two respects. First, the additivity of 
the counting measure allows one to appeal to an ergodic theorem, whereas subadditivity requires 
the use of subadditive ergodic theorems (such that of Akcoglu and Krengel [AK81]), and the 
ergodicity of the point process itself (which is used in a weaker form in [PenOl]). Second, and 
more importantly, the additivity of D i-> ((D) and the uniform bound < ((D) < C\D\ 
imply that the contribution of any subset D' of D to S(D,() is uniformly bounded by C\D'\ 
— which is crucial in [PenOl]. In the subadditive case, this does not hold: The contribution of 
a subset D' of D to S(D, () is not a priori bounded by C|D'|. This compels us to appeal to 
the stabilization properties introduced by Schreiber, Yukich, and the second author in [SPY] — 
which are complementary to the ones in [PenOl]. The rest of the proof relies on a percolation 
argument and the Borel-Cantelli lemma. 

In the second part of this work we apply this general result to classical Euclidean optimiza- 
tion problems on the random parking measure. In particular, we shall prove a so-called umbrella 
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theorem, which allows to cover at once a wide class of problems, including the traveling sales- 
man, the minimum spanning tree, and the minimal matching. This is one of the first examples 
where subadditive arguments are combined with stabilization properties. We refer to the recent 
survey of Yukich [YuklO], where both types of arguments are presented and used — although 
not combined. 

The third part of this work is concerned with another application of our general result: the dis- 
crete model for rubber introduced in [GLTV], whose thermodynamic limit is studied in [ACG]. 
In particular, as recalled in Subsection 4.1, the model under investigation is based on the notion 
of stochastic lattices (in the sense of random point processes) and on the associated Delaunay 
tessellations. In a nutshell, every edge of the tessellation represents a polymer chain (the vertices 
are then permanent cross-links). As argued in [GLTV] and [ACG], it is reasonable at first order 
to consider cross-links at zero temperature and the polymer chains at finite temperature. The 
free energy of the E-rescaled network of polymer chains in a domain D of IR 3 is then given by 
the sum of the free energies of the edges, plus a volumetric term accounting for the incompress- 
ibility of the network (the volume of each simplex is conserved). Under some assumptions on 
the energy terms (which are standard in the theory of homogenization of integral functionals), 
and under assumptions on the stochastic lattice compatible with the random parking measure, 
it is proved in [ACG] that the free energy functional of the e-rescaled network in D (seen as 
a function of the position of the cross-links) T-converges (as e goes to zero) to a continuous 
energy functional of the type 



«>->■/ H / hom(Vu(a;))dx, 
Jd 

where the "homogenized energy density" Whom is a quasiconvex frame-invariant function, as 
encountered in continuum mechanics. The map u : D — > M 3 is a deformation. In particular, 
this type of convergence ensures that minimizers of the free energy of the discrete system con- 
verge (up to extraction) to minimizers of the continuous energy functional — which yields a 
rigorous derivation of nonlinear elasticity compatible with minimization. In addition, if the sto- 
chastic lattice is statistically isotropic, then the associated homogenized energy density Whom 
is isotropic, as expected in rubber elasticity. The stochastic lattices considered in [ACG] satisfy 
the following three properties: ergodicity, non-empty space condition, and hard-core condition. 
Two questions were left open in [ACG]: 

• Do there exist such stochastic lattices which are indeed statistically isotropic ? 

• What happens to the thermodynamic limit in D if the stochastic lattice on R d is replaced 
by an approximation on D (in which case boundary effects may appear and rule out 
stationarity) ? 

This article gives a clear answer to both questions. As we quickly show in Subsection 2.1, the 
random parking measure in M. d is an example of stochastic lattice which is statistically isotropic, 
ergodic, and satisfies the non-empty and hard-core conditions. In Subsection 4.2, we shall prove 
that the model studied in [ACG] can be recast in terms of a subadditive function satisfying the 
assumptions introduced above, so that the identity corresponding to (2) will hold true. This 
allows us to complete the program initiated in [ACG] both in terms of lattices generated in D 
instead of M. d (which seems more reasonable from a physical point of view), and in terms of 
numerical approximations (effective computations for Whom i n [GLTV] are based on random 
parking in bounded domains Qr). 

This paper is organized as follows. Section 2 is dedicated to the study of the qualitative 
properties of the random parking measure, and to the proof of (2). In Section 3 we apply 
the main result to classical Euclidean optimization problems. Finally, Section 4 is concerned 
with the application of the main result to the discrete model for rubber studied in [ACG]. In 
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particular, we show that this model satisfies the assumptions of the main result, which completes 
the analysis of the model when the stochastic lattice is the random parking measure. 

We make use of the following notation throughout the article: 

• d, n > 1 denote dimensions; 

• M+ = [0,+oo); 

• 0(R d ) is the set of bounded nonempty Lipschitz subsets of R d (a bounded measurable 
set D, whose boundary is denoted by dD, is Lipschitz if there exists a finite collection 
of relatively open sets [/& C dD such that dD = U^U^, and Uk is the graph of a 
Lipschitz function on R d_1 up to a rotation); 

• For all D G 0(R d ), \D\ denotes its d-dimensional Lebesgue measure, dD denotes its 
boundary, and D its closure; 

• For all D G 0(R d ) we set 0{D) = {A G 0{R d )\A C D} (note that D G 0(D)); 

• For all x G R d , |x| is the Euclidean norm of x, and cfe denotes the Euclidean distance, 
i.e. d,2{x, A) := inf{|x — y\ : y G A} for nonempty A C R d . 

• diam denotes the Euclidean diameter of a subset of R d ; 

• M. nxd is the space of n x d matrices, that we simply denote by M. d when n = d; 

• For all A G M nxd , we define <p\ : R d ->• R n by x i-> Ax, and denote by |A| the 
Frobenius norm y / Trace(A T : A) (which is the operator norm of (p\ associated with 
the Euclidean norms of R d and R n ); 

• SOd is the set of rotations in R d ; 

• Ai and Afr denote the spaces of finite and locally finite point sets in R d , respectively; 

• For all open D C R d , C(D,R n ) and C$°(D,R n ) respectively denote the space of 
continuous functions from D into lR n , and the space of smooth functions which have 
support in D \ 3D; 

• For all D G 0(R d ) open, n > 1, and p G [l,oo], LP(D,R n ), W^ p (D,R n ), and 
W ' P (D, R n ) denote the Lebesgue space of p-integrable (or essentially bounded if p = 
oo) functions from D into R n , the Sobolev space of p-integrable functions on D whose 
distributional derivatives are p-integrable (or essentially bounded if p = oo), and the 
closure of C^(D, R n ) in W^ P (D, R n ), respectively; 

• When no confusion occurs we use the short-hand notation C(D), L P (D) etc. for 
C(D,R n ) in LP(D, R n ) etc.; 

• For all D G 0(R d ) open and u G W 1,oc (D), we denote by ||«||Lip the Lipschitz 
constant of u on D. 

2. Random parking and subadditive ergodic theorems 

2.1. Random parking. We first quickly recall the graphical construction of the random park- 
ing measure £ in a Borel set A C R d . Let V be a homogeneous Poisson process of unit 
intensity in R d x R + . An oriented graph is a special kind of directed graph in which there is no 
pair of vertices {x, y} for which both (x, y) and (y, x) are included as directed edges. We shall 
say that x is a. parent of y and y is an offspring of x if there is an oriented edge from x to y. By 
a root of an oriented graph we mean a vertex with no parent. 

The graphical construction goes as follows. Let po > and let B denote the Euclidean ball 
in R d of radius po centred at the origin. Make the points of the Poisson process V on R d x R + 
into the vertices of an infinite oriented graph, denoted by Q, by putting in an oriented edge 
{X, T) ->• (X', T') whenever (X' + B) n (X + B) / and T < T'. For completeness we also 
put an edge (X, T) ->• (X 1 , T') whenever (X' + B) n (X + B) ^ 0, T = T', and X precedes 
X' in the lexicographical order — although in practice the probability that V generates such 
an edge is zero. It can be useful to think of the oriented graph as representing the spread of 
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an "epidemic" through space over time; each time an individual is "born" at a Poisson point in 
space-time, it becomes (and stays) infected if there is an earlier infected point nearby in space 
(in the sense that the translates of B centred at the two points overlap). This graph determines 
which items have to be accepted. 

For (X, T) G V, let C( X ,T) (the "cluster at (X, T)") be the (random) set of ancestors of 
(X, T), that is, the set of (Y, U) G V such that there is an oriented path in Q from (Y, U) to 
(X, T). As shown in [PenOl, Corollary 3.1], the "cluster" C {X ,T) is finite for (X, T) e V with 
probability 1. It represents the set of all items that can potentially affect the acceptance status of 
the incoming particle represented by the Poisson point (X, T). The set of accepted items may 
be reconstructed from the graph Q, as follows. Let A C M. d be a (possibly unbounded) Borel 
set, and let Va denote the setVH (A x M + ), i. e. the set of Poisson points that lie in A x R + . Let 
Q\A denote the restriction of Q to the vertex set Va- Recursively define subsets Fi(A), Gi(A), 
Hi(A) of A, % = 1, 2, 3, . . . as follows. Let F\ {A) be the set of roots of the oriented graph Q\ A , 
and let Gi(A) be the set of offspring of roots. Set Hi(A) = F 1 (A) U G 1 (A). For the next 
step, remove the set H\(A) from the vertex set, and define F2(A) and G 2 (A) the same way; so 
F 2 (A) is the set of roots of the restriction of Q to vertices in Va \ H\{A), and G 2 (A) is the set 
of vertices in Va \ #1 (A) which are offspring of F 2 (A). Set H 2 (A) = F 2 (A) U G 2 (A), remove 
the set H 2 (A) from Va \ H\{A), and repeat the process to obtain F 3 (A), G 3 (A), H 3 (A). 
Continuing ad infinimum gives us subsets Fi(A), Gi(A) of Va defined for i = 1, 2, 3, ... . 
These sets are disjoint by construction. In the case when A = ~R. d , we drop the reference to A 
and use the abbreviation F { and G t for Fi(R d ) and Gi(R d ). 

As proved in [PenOl, Lemma 3.2], for every bounded nonnull Borel set A in R d , the sets 
Fi(A), Gi(A), F 2 {A), G 2 (A), . . . form a partition of V A , and F lt G x , F 2 , G 2 ,... form a 
partition of V with probability 1. 

Definition 2.1. The random parking measure £ in A is given by the projection of the union 
U ( j*i 1 Fi(A) on M. d . Likewise, the random parking measure £ in M. d is given by the projection of 
the union Ug^F on R d . 

Definition 2.2. Let S C R d be a locally finite set of points. We say that £ is general if no d + 1 
points lie in the same hyperplane and if no d + 2 points lie in the same hypersphere. 

Definition 2.3. Let < p\ < p 2 < 00. Given D C M. d , suppose E is a subset of D. Then 
S is said to be (pi, p 2 )-admissible in D iff for all x 7^ y G S, \x — y\ > p\, and for all 
z G D, B{z, p 2 ) flS/8, where B{z, p 2 ) denotes the open ball centred in z of radius p 2 (and 
B(z, 00) := R d )). Let A pltP2 denote the class of general point sets that are (pi, p 2 ) -admissible 

inIR d . 

In other words, an admissible point set in D is one which satisfies the hard-core and empty 
space conditions. We shall sometimes write simply 'admissible' for 'admissible in M. d \ 

Definition 2.4. A random subset £ of M d is called a point process. We say that a point process 
£ in M. d is isotropic if the distribution of £ and of R£ are the same for every rotation R G SOd- 

We are now in position to prove that the random parking measure £ in U. d (which is admissible 
by construction) is an isotropic admissible stochastic lattice in the sense of [ACG]. 

Proposition 2.1. The random measure £ is stationary (under real shifts), ergodic, isotropic, 
and almost surely general. 

Proof of Proposition 2.1. We prove each property separately. 

Step 1. Stationarity. 

By definition, the Poisson point process V on M. d x M + and its translation (xi, . . . , xj,, 0) + V 
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have the same distribution for all x = (xi, . . . , Xd) G R d . Hence the graphical construction is 
stationary, and therefore also the random parking measure £. 

Step 2. Isotropy. 

The proof of the isotropy of £ is similar to the proof of the stationarity. For all R G SOd, the 
Poisson point process V and its rotated version HV := {(Rx, t) : x G R d , t G M + , (x, i) € 
V} have the same distribution, which implies that the random parking measure £ is isotropic. 

Step 3. General position. 

For all t > we set P t = {x G R d : 3r € [0, t], (x, r) G £>}. Since £ C U neN P n , and since 
V n C 7-Vn for all n G N, the event that £ is not in general position is contained in the union 
over n G N of the events that V n is not in general position. Since V n is a Poisson process of 
intensity n in R d , the probability that V n is not in general position is zero, so that the union of 
this countable set of events has also probability zero, and £ is almost surely general. 

Step 4. Ergodicity. 

Let us view £ and V as elements of Af f and A^ 1 respectively. Also let us extend V to a 
homogeneous Poisson process of unit intensity on the whole of M. d+1 (also denoted V). Let T x 
denote translation by an element x of M. d (acting either on A d ? or A\ ^ according to context). 
Then as described earlier in this section, £ is the image of V under a certain mapping h from 
A d i l to Afp which commutes with T x for any x G R d , that is T x o h = h o T x . 

Suppose A is a measurable subset of Aft which is shift-invariant, meaning T X (A) = A for 
all x G R d . Then for x G R d , 

h-\A) = h-\T x {A)) = T x {h~ 1 A) 

so h^ 1 (A) is invariant under the mapping T x acting on A^ 1 . Now, V is ergodic under 
translations, that is if B C Aft satisfies T X (B) = B for some non-zero x G M d+1 , then 
P[B] G {0, 1}. See for example the proof of Proposition 2.6 of [MR]. Therefore with A as 
above, 

P[£ eA] = P[V G h- l {A)} G {0, 1}. 

Thus £ is ergodic. □ 

2.2. Main result. For i G Z d , let C { := i + [-1/2, l/2] d be the closed rectilinear unit cube 
centred at i. 

Let < pi < p 2 . Given D G 0(R d ), and R > 0, let D R := {Rx, x G D} be the dilation 
of D by a factor R > 0. For r > set D R)T := {x G D R : d 2 (x,IR d \ D R ) > r). Also set 
D Rfi := D R , and let V(D) := {D R:r : R > 0, r > 0, D^, G C(M d )}. The next result shows 
that D Rir G T>{D) whenever r/R is small enough. 

Lemma 2.2. Suppose D G C(R d ). 77je« ?/jere exwto 5 > .rac/i f/iaf L> M G C(IR d )/or a// 
iG]0,<5[. 

Proof. Assume without loss of generality that d > 2 and D is closed. For r > let i^, denote 
the closed Euclidean ball of radius r centred at the origin in M d_1 . Suppose x G dD. By the 
definition of 0(R d ) we may assume after translation and rotation that x = and that there exist 
r > and a > 0, and Lipschitz function g : B' 2r — > R such that |g(u)| < a for all u G B' 2r and 

Dn{B' 2r x [-4a, 4a]) = {(u, y)|u G B' 2r ,y G [-4a, 5 (u)]}. 

Suppose < t < min(r, a). Then 

D lyt C\(B' r x [-2a, 2a]) = {(u,y)\u G B' r ,-2a<y < inf 5 w (n)} 



RANDOM PARKING, EUCLIDEAN FUNCTIONALS, AND RUBBER ELASTICITY 



where for w G B' t we set g w {u) = g{u + w) — y/t 2 — \w\ 2 . Since the infimum of Lipschitz 
functions with common Lipschitz constant is itself Lipschitz, it follows that the intersection of 
dDij with the interior of B' r x [—2a, 2a] is the graph of a Lipschitz function of d — 1 variables. 
We may then use a compactness argument to deduce that there exists 5 > such that for 
< t < e the set D\t is Lipschitz. □ 



It follows from Lemma 2.2 that given any D G 0(R ) and < a < 1, the set Dr^rc is 
in 0(M. d ) for all large enough R. We now define some properties of functions parametrized by 
point sets (or restrictions of point sets on bounded domains). 

Definition 2.5. Let D G 0(R d ). A measurable function S : V(D) x A PltP2 -> R is said to be: 

• local on V{D) if for all D G V(D), and all C, C € A PuP2 such that C n D = ( D A we 
have 

s(AO = s(AO; 

• insensitive to boundary effects on T>(D) if there exists < a < 1 such that we have 
(3) hmsup sup < —j ; > = 



• Also, S 1 is said to have the averaging property on A pijP2 with respect to D if for any sta- 
tionary and ergodic random lattice ( whose realization almost surely belongs to A pitP2 , 
there exists a deterministic S = 5(C) G R such that almost surely 

(4) lim ^^1 = S. 

R^+oo \D R \ 

If S is local, then we can properly define S(D, C) for a point set ( only defined in D G 
0(M. d ), provided it is general and (pi, ^-admissible in D. We do this as follows: let (* be a 
choice of general and (pi, ^2) -admissible set in M rf such that ("* H -D = C> an d set S(D, Q := 
S(D, (*). By locality this definition does not depend on the choice of (*, and there does exist at 
least one such choice of C* because if we choose po £ (pi , P2) and perform the random parking 
process in R. d as described in section 2.1, but rejecting any point that lies within distance po of 
any point of f , then the resulting set of accepted points, together with the point set ( itself, will 
be (pi, P2) -admissible on the whole of R , and almost surely general. 

In Sections 3 and 4 we shall give two examples — or rather two classes of functionals — 
which satisfy the above properties. Both examples exhibit subadditivity, and we shall make 
use of the Akcoglu-Krengel subadditive ergodic theorem to show they satisfy the averaging 
property. 

Let < pi < po < p2 < 00, and for D G 0(M. d ), let £ (respectively £ D ) be the random 
parking measure on M. d (respectively on D) with parameter po (as in Definition 2.1). The main 
result of this section is the following: 

Theorem 2.3. Let D G 0(M d ). If the measurable function S : T>(D) x A P1 , P2 — > M + is local 
on T>{D), insensitive to boundary effects on T>{D), and has the averaging property on A pitP2 
with respect to D, then with S given by (4), almost surely 

^ ,. S(Dr^ Dr ) ,. S(D R ,Q -= 

(5) hm — — | = km — — — — = S. 

R^+oo \Dr\ _R->+oo \Dr\ 

Theorem 2.3 can be seen as a version of the 'subadditive ergodic theorem' for random park- 
ing measures on homothetic sets. It is proved using the following result due to Schreiber, Yukich 
and the second author [SPY], which gives an 'exponential stabilization' property of £. Let us 

say y C R d x R+ is temporally locally finite if y n (R d x [0, t\) is finite for all t G (0, 00). 
Also, for r > let B r denote the Euclidean ball {2; G l rf : |x| < r}. 
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Lemma 2.4. There exists a nonnegative random variable V such that (i) V has an exponen- 
tially decaying tail, i.e. for some constant c > the probability that V exceeds t is at most 
cexp(— c~ l t) for all t > 0; and (ii) for all temporally locally finite y C (M d \ By) x M + 
the measures £ Bv and ^ Bv ^ coincide on [—1/2, l/2] d , where ^ Bv ^ is the parking measure 
induced by the union ofPf] (By X M + ) and y. 

Proof. Essentially this result is Lemma 3.5 of [SPY]. That result says only that the total number 
of points in [—1/2, l/2] rf coincide for ^ B ^y and ^ Bv , not that the measures themselves coincide 
on [—1/2, l/2] rf . However, the proof in [SPY] in fact demonstrates the stronger statement given 
here. The proof in [SPY] assumes that po < 1/4 (see [SPY], p. 173) but once we have the result 
for this case, we can deduce it for the general case by a scaling argument. □ 

In the next lemma, we use Lemma 2.4 to show that S(D RiRa , £) is a good approximation to 
S(D R>Ra ,^). 

Lemma 2.5. Suppose D G 0(R d ) and S : V(D) x A PltP2 ->■ M + is local on V(D). Suppose 
< a < 1. Then there exists an almost surely finite random variable Rq such that for R> Rq 
we have S(D R , Ra ,£ Dr ) = S(D R , Ra , £). 

Proof. For i G 7L d and r > 0, let B\ denote the closed Euclidean ball of radius r centred 
at i, and let Qi/ 2 denote the rectilinear unit cube centred at i. By Lemma 2.4, there exists a 
family of random variables (V(i),i G Z d ) which are identically distributed with exponentially 

B^ V B^ — i 

decaying tails, such that £ v w coincides with £ V W on Q\i 2 , for all temporally locally finite 
y C (R d \ By,*}) x R+, where £ v (»)' denotes the parking measure induced by the union of 

V n (B v(i) x R+) and y. 

Let D RRa := {i G Z d : Q\ /2 H £>#,#« / 0}. For fc G N, let E(k) denote the event that 
there exists R G [fc, fe + 1[ and i G D R Ra such that By« \ D R ^ 0. If E{k) does not occur, 
then for all R G [A;, k + 1[ and all i G DJj fia we have By,* C -Dr, and therefore both £ and 

£ R coincide with £ v w on Q\i 2 . Hence £ and £ R coincide on the whole of D R<R a. 

Choose K such that D C B K . Suppose k G N and i? G [k,k + 1[. Then for i G D RRa , 
since i is distant at least R a — d from M. d \ D R by the definition of D R Ra , we have that 
B\ a _ d C B l Ra _ d C Dr. Also, we have D RRa C Bjnj+d C B K{k+1)+d . 

Therefore, for large enough k G N the probability that E(k) occurs is bounded by the cardi- 
nality of 7L d fl B K ^ k+ ij +d , multiplied by the probability that V(0) > k a — d, and by Lemma 
2.4 this is bounded by ck d exp(— c~ l k a ) for some universal constant c > 0. Hence, the prob- 
ability of E{k) is summable over k G N, and from the Borel-Cantelli lemma we deduce that 
almost surely, E(k) occurs for only finitely many k so that there exists Rq such that for all 
R > Rq the measures £ and ^ Dr coincide on D RjR a, so that S(D RtR a,£) = S(D RtR a,£ DR ) 
by locality. □ 

Proof of Theorem 2.3. Since £ is ergodic by Proposition 2.1, and S has the averaging property 
in the sense of Definition 2.5, there exists S G M + such that (4) holds, i.e. the second equality 
of (5) holds. Choose a G (0, 1) such that (3) holds. To prove the first equality of (5), for all 
R> 1, we write £ R for £ Dr and rewrite the averaged energy as 

,~ S(Dr,Z r ) S(D R ,Q 

(6) — ^n — = |n ; + ei(-R) + e 2 (R) + e 3 (R), 
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with 



ei(R) 
e 2 (R) 
e 3 (R) 



S(Dr,r«,0 S(D R ,Q 
\Dr\ \D R \ ' 

S(Dr,r«,Z R ) S(D R , Ra ,Q 
\Dr\ \D R \ ' 

S(D R ,£ R ) S(D R>Ra ,C R ) 



\Dr\ \Dr\ 

Now e 2 (R) tends to zero almost surely by Lemma 2.5, and by the insensitivity to boundary 
effects e\(R) and e%{R) tend to zero almost surely as R — > oo, so that the first equality of (5) 
follows. □ 

3. Application to classical Euclidean optimization problems 



A measurable functional S : 0(Mr) x Aff — > R is said to be translation invariant if for 
every y € R d , £> € 0(R d ), and every ( € Af f , we have S(D + yX + y) = S(D, (), where 
for any subset U of R d , U + y := {z + y : z£ U}. 

Definition 3.1. Let S : 0(R d ) x Aff — > R be a translation invariant measurable functional, 
and let p > 0. We say that 5 is 

" C ' V ~ ' >ULH llltU ,U1 UJJ L 3 --Y/- 



localized if there exists 5 : Af -> R such that for all ( G .4f f and D e O 
S(D,0 = S((nD), 



• almost subadditive of order p if there exists (7 > such that for all Ci, C2 £ Aft and all 



D £ 

S(ACiUC 2 ) < 5(ACi)+5(ACj)+Cdiam(I>f; 

• (strongly) superadditive on rectangles if for any partition of a semi-open rectangle D 
into a finite number k of semi-open rectangles {-Di}ie{l,...,fc}> 

fc 

s(ao > X>(a,o; 

• smooth of order p if there exists C > such that for every rectilinear cube D and all 

(1,(2 eAf f , 

(7) |5(AClUC2) " S(DXi)\ < Cdiam(D) p (card( 2 nD) ( - d - p V d . 

• homogeneous of order p if for all a > 0, D € C(R rf ), and £ G Af/-, we have 
S(D a , a() = a p S(D, (), where a( := {ax : x e Q. 



Definition 3.2. Two measurable functionals Si and S 2 are said to be pointwise close of order 

id 
K lf> 



p > if for every rectilinear cube D C M. d and every ( € «4f f , we have 



1^(0,0-52(0,01 / ,,„„,(,(-,)/,, 



<8) diam(oy '- =H' Card ^ P » 

in the sense that there is a function g : Z + — > R + with lim^oo k( p ~ d ^ d g(k) = 0, such that the 
left side of (8) is bounded by g (card (nD) uniformly over D and £. 

Regarding terminology, note that the condition that S be localized (in the sense of Definition 
3.1) is stronger than the condition that S be local (in the sense of Definition 2.5). 

A translation invariant measurable function S that is also homogeneous of order p > and 
satisfies S(D, 0) = for all D, is called a Euclidean functional. See [Yuk98] for examples of 
Euclidean functionals. 
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We are now in position to state an umbrella theorem for translation invariant functionals on 
the random parking measure. Let £ and ^ D be as in Definition 2. 1 

Theorem 3.1. Let po > 0. Let 1 < p < d, and let S and T be measurable translation invariant 
smooth functionals of order p. Assume that S is localized and almost subadditive of order p, 
and T is superadditive. If S and T are pointwise close of order p, then there exists deterministic 
S £ M + such that for all D £ Q(R ), we have almost surely 

lim S(D R ,S°*) = Hm S(D R ,Q = -g 

iJ^oo \D R \ R-+00 \D R \ 

Remark 1. We do not assume homogeneity of S and T in the statement of this theorem. How- 
ever, if in fact we do assume that S and T are homogeneous of order p, then to check S is smooth 
of order p we need to check (7) only for D = [0, l] d , and to check S and T are pointwise close 
we need to check (8) only for D = [0, l] d . 

Remark 2. By [Yuk98, Lemmas 3.2, 3.5, and 3.7], the p-power weighted total length version 
of the traveling salesman problem, the minimal spanning tree, and the minimal matching all 
satisfy the assumptions of Theorem 3.1 with exponent p. In particular, these examples are ho- 
mogeneous and the associated superadditive Euclidean functions T are the so-called canonical 
boundary functionals (denoted by S R in [Yuk98]). 

We believe there could be interesting examples where homogeneity fails but our theorem is 
still applicable, for example to do with the number of percolation components. See [Yuk98], 
page 47 for a discussion of examples where homogeneity fails. 

Remark 3. The proof of Theorem 3. 1 shows that its conclusion holds for any bounded connected 
subset D of R^ whose boundary has zero (i-Lebesgue measure. 

Remark 4. In the proof of Theorem 3.1 we show that S has the averaging property on A pij00 
with respect to D for any D £ 0(M. d ) and for any p\ > (in particular we need only the 
hard-core property). 

Proof of Theorem 3.1. We split the proof into five steps. We first prove that S satisfies the 
insensitivity to boundary effects. In the second step, we address the averaging property of T on 
cubes with vertices in Z d , then that of S on cubes with vertices in M. d , and finally on any domain 
in 0(K d ). Appealing then to Theorem 2.3 we will conclude the proof. 

Step 1. Insensitivity of S to boundary effects. 

Let D £ 0(R. d ) and < a < 1. Let Q be a cube aligned with the canonical basis of M. d and 

which contains D. By localization and smoothness of S, for all £ £ A Pl! oo, 

\S(D R ,0-S(D R , Ra ,0\ = \S(Q R ,(nD R )-S(Q R ,(nD RRa )\ 

< CR p ( ca vdCnD R \D R ^ d -^ d . 

Since the boundary of D has zero (i-Lebesgue measure, we can cover the boundary of D R by 
^(i-a)d^ balls of radius R a . Provided R > 1, the corresponding balls of radius 3R a cover 
D R \D RjR a and by the hard-core condition, each such ball contains 0(R da ) points of (. Hence 
cardC n D R \ D RjRa is o(R d ), and hence \S(D R , ()-S(D RRa ,()\ = o{R d ). Thus (3) holds, 
so S is insensitive to to boundary effects on V{D). 

Step 2. Superadditive ergodic theorem on rectangles with vertices in Z d . 

Let C be a stationary ergodic random lattice taking values in A Plt00 - Let 1Z denote the collection 

of half-open rectilinear rectangles in M. d with vertices in Z d , and let C be the class of all cubes 
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in 1Z. By localization S(D, 0) is a constant 5(0) independent of D. By pointwise closeness, 
there is a constant C such that for every DeCwe have that 

T(D, 0) < S(D, 0) + C(diam(D)) p 

and hence by taking £i = and £2 = C i n me definition of the smoothness of T, we have for a 
possibly different C that for all D G C, almost surely 



T(D,Q < Cdiam(L>) p (l + cardCnL») 



(d- P )/d 



The hard-core condition combined with an elementary geometric argument imply there exists 
C > depending only on p\ such that card £ n D < C\D\ for every D £7Z. Hence, there is a 
further C such that if D G C then 

(9) ET(AC) ^ Cl^l- 

However, this implies that (9) also holds whenever D £ 1Z, since if not, then we could take a 
large cube in C which was the union of disjoint translates of D and using stationarity of ( and 
translation-invariance and superadditivity of T, we would have a violation of (9) for this cube. 
We may now apply the Akcoglu-Krengel superadditive ergodic theorem [AK81] to the func- 
tional D \— > T(D, £) which is stationary under integer shifts. By the Akcoglu-Krengel theorem, 
for any unit cube Q G C there is a random variable S(Q) such that almost surely 

(10) lim ^%i^ = 5(Q). 

R->+oo,RGN R d 

Indeed, if Q is contained in the positive orthant then the sequence of cubes Qr is regular in the 
sense of Akcoglu and Krengel. If not, then (since it is a unit cube) Q is contained in some other 
orthant and the argument is similar. 

By the pointwise closeness of 5 and T, and the hard-core condition on £, (10) implies that 
f or Q G C a unit cube we have almost surely 

(ID lim ^%^ = 5(Q). 

We claim that S(Q) is deterministic and the same for all unit cubes Q G C (hence subsequently 
denoted by simply 5). Indeed, given such a Q and given z G W 1 , we may deduce by localization 
and smoothness of 5 that 

\S(Q R ,Q-S(Q R n(Q R + z),0\ < CRP(cardCDQ R \(Q R + z)f' p)/d . 

By the hard-core property of £, there is a constant C depending on p\ and z such that 
card C H Q R \ (Q R + z) is asymptotically bounded by CR d ~ l , so that 

n ~ r \S(Q R ,C)-S(Q R n(Q R + z)X)\ n 
(12) lim sup — -, = 

R ^oo,nen R 

and likewise R" d \S(Q R + z, () - S(Q R n (Q R + z), ()\ -> as # -> 00 through N. Hence 
by (11) we have 

lim *(3* + *>fl = 5(Q). 

Hence S(Q) is shift invariant, and therefore constant by ergodicity Also, by translation invari- 
ance of 5 and stationarity of (, the distribution of S(Q) is the same for all unit cubes Q G C, so 
we have justified our claim. 
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Step 3. Extension of (1 1) to general cubes. 

Let C be as in Step 2. Let Q G C be a cube of side k with k € N. Let Q 4 , 1 < i < k d , be disjoint 

unit cubes in C, whose union is Q. Then by superadditivity of T, 

fed 

(13) liminf r( ^ ,C) > liminf ^Wlfi.C)^^ 

i=i 
Also, by repeated use of the almost subadditivity and localization of 5, as in [Yuk98, Remark 1 
(3.5)], there is a constant C(k) such that 

k d 



S(Qr,C) < J>((Q J ) fi ,C) + CW, 



so that 

fed 

(14) limsup S{QR d ' C) < limsup -L £ S((Q%, C) = k d S. 

2 = 1 

Moreover by closeness of T and S, and the hard-core property of £, 

(15) limsup ij|T(g fl ,C)-S , (Qfl,C)l = 0. 
R-+oo,R€H n 

Combining (13), (14) and (15) we may deduce that (10) and (11) hold for all Q £ C, with 

5(Q) = \Q\S. 

Next, for Q € C we relax the assumption that R € N in the asymptotic formula (11). 
Denoting by [R] the integer part of R, by using the localization and smoothness of 5 and the 
hard-core property of ( as in the proof of (12) above, we have 

\S(Qr,Q-S(Q [r] ,C)\ n 
lim sup — 3 = U. 

Combined with (11) this shows that 

(16) lim ^9jh0 = s. 

K^+oo \Q R \ 

Now let Q be a general rectilinear cube in R. d . For all R > large enough, there exist k R € N 
with k R > R — 2 > and a rectilinear cube Q k R with vertices in 7L d and sidelength k R such 
that Q fc fl C Q R . As above, the hard-core property of £ yields the estimate 

card C n Q RR \ (Q k A) R = card C n (Q R \ Q k *) R < CR d R d - 1 
so that localization and smoothness of S yields 

limsnp |s((^) R ,Q- 5(^,01 £ oJ8M/<i 

so using (16) applied to Q k R, and the arbitrariness of R, we obtain 

(17) lim ^Jh0 = s. 

R-++oc \Q R \ 

Step 4. Extension of (17) to general domains D G 0(R. d ). 

Let £ be as in Step 2. Let D + be a cube aligned with the canonical basis of Z d , containing D, 

and such that its sidelength 5 is bounded by diam(D). 
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We proceed by approximation. For every integer j > 1, let Ej = 2~ 3 5, and split D + into 2° 
cubes {Q J ' J }ie{i ... 2-? 1 *} °f sidelength Sj. Define D 3 = \J ieI jQ : >'' 1 as the union of those cubes 
Q 3 ' 1 (indexed by i £ P) whose intersection with D is not empty, and define D 3 = U ie jjQ 3 '' 1 as 

the union of those cubes QJ> i (indexed by i £ P) which are contained in D. In particular, P is 
finite and D j CDc D 3 C D+ . 

Let rj > 0. Using the fact that the boundary of D has zero Lebesgue measure, choose j such 
that \D 3 \ D 3 | < rj, so that D 3 \ D 3 consists of at most ijej cubes of side Sj. This j will be 
fixed for a while. By the hard-core constraint there is a constant C depending only on d and p\, 
such that for R so large that Rej > 1 we have 

(18) card C n {D 3 \ D j ) R < Cr ] e^ d (Re j ) d = CR d rj. 
Hence by localization and smoothness of 5, 

y \S(D R ,C)-S((D 3 ) R ,()\ |5(D+,cng fl )-g(D+ cn(jJJ) fl )| 

hm sup — d = ton sup * —^ 

(19) < Cr]( d - p V d . 

Using localization and repeatedly using almost-subadditivity, as in [Yuk98, Remark 1 (3.5)], 
yields for some constant C(j) that 



,i 



V 

S((d 3 ) r ,() = S(D+,cn(D*) R ) < J2s((Q j nR,(n(D 3 ) R ) + c(j)R p - 

2=1 

Since 5((Q-?> i ) ii , ( n (D 3 ) R ) is a constant 5(0) for i <£ P, this turns into 
S((D 3 ) R ,() < j;#%C) + C(j')F. 

Hence, using (19) and using (17) for each cube Q 3,i , i £ P, this yields 

(20) ton sup R~ d S(D R X) < Crf^r + ^ \Q j ' { \S < Crf^r + (|D[ + rj)S. 
R->+oc . £7i 

We now turn to the lower bound. By almost subadditivity and localization of 5, 

(21) S(D+,() < S(D R , C) + S(D+ \ D R , C) + CRP. 
By localization and smoothness of 5, 

\S(D+\D R ,Q-S(D+\(D?) R ,Q\ < CR?(caxd(nD R \(D%)( d -PV d , 
and by (18) this gives us 

(22, Vmmp VHD&\I>*,0-SlDi\<B»)*,a\ < c „^. 

R->oa R 

We now treat S(D^ \ (D?)r, C) : by localization and almost-subadditivity, 
S(D+\(D?) R ,0 = S(D+,(nD+\(D 3 ) R ) 

2> d 



< ^5((Q^,cn J D+\(^) fl ) + C(i)^ 

2=1 

(23) = J^S((Q 3 >%,C) + C(j)R P - 



HP 
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Applying (21), (22) and (23) in turn, and using (17) for each cube Q 3 ' 1 and for D + , we obtain 
liminf R~ d S(D R ,() > \D + \~S - limsup R- d S(D+ \ D R ,() 

> \D + \S - CT] l - p/d - lim sup R- d S{D+ \ {D j ) R , () 

> \D + \S-CT] l - p/d -\D + \D j \S 

> (\D\- v )s-Cn l - p/d . 

Combined with (20) and using the arbitrariness of r/, this concludes the proof of the averaging 
property (4) for S. 

Step 5. Conclusion. 

We can now apply Theorem 2.3 to S, since in Step 1 we have shown the insensitivity to bound- 
ary effects, and in Step 4 we have shown the averaging property. □ 

4. Application to the derivation of rubber elasticity 

4.1. A discrete model for rubber and its thermodynamic limit. Let < p\ < p2 < oo. 

Suppose that £ E A PljO0 and C has at least d + 1 elements (so its convex hull has strictly 
positive (i-Lebesgue measure, since C is assumed general). Then there is a unique Delaunay 
triangulation of the convex hull of £, by simplices with edges given by the edges of the Delau- 
nay graph of C, which is itself the dual graph of the Voronoi tessellation (see [OBSNC], and 
[Del34] for Delaunay tessellations of R d , and for instance [For95] for Delaunay tessellations of 
a bounded domain). By convention we consider simplices to be open sets. Let T{C) denote the 
Delaunay triangulation of the convex hull of C. For any d+ 1 points of C, the simplex generated 
by these points is in T{C) if and only if no point of £ lies in the interior of the circumsphere 
of these points. We denote by M(C) the associated neighbour pairs, that is, those unordered 
pairs of points {x, y} such that (x, y) is an edge of T := T{C). For all e > and all open sets 
D € 0(M. d ), we define a space of continuous piecewise-affine functions &®(£) on D, by 

(24) ©f (£) := {u E C{D,W n ) \ \/T G T(£), witheT n D ± 0, U| eTnD is affine}. 

Let D denote the closure of D. From now on, we identify u : eC n D — > W 1 with its class of 
piecewise-affine interpolations (still denoted by u) in &®(£) C W 1 '°°(D,M. n ). Note that the 
extension of u : eC n D — > W 1 to D \ UxeT,TcDT is not uniquely defined — as we shall see, 
the energy under consideration does not depend on the extension. In order to define an energy 
functional on the set &^(£), we first introduce the following energy functions: 

Definition 4.1. Let p > 1. We denote by U p the subset of functions f nn of C(R d x R n , R+) 
for which there exists C > such that, for all z G R d and s G R n , 

(25) ^\s\ p -C<f nn (z,s)<C(\s\ p + l). 

We denote by V p the subset of functions W vo \ of C(M nxd , R + ) for which there exists C > 
such that for all A G M nxd , 

(26) W vol (A) < C(\A\P + 1). 

Let p > 1 and f nn G Up, W vo \ G V p . For all u G L X {D, R n ) and open D G 0{R d ) we then 
set 

(27) F D (C u) = { F nn,e{£,u) + Fg l£ {C,u) if U G 6? (£), 

£ y ' ' ' I +00 else, 
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where we set 

(28) < £ (/» = £ e d f nn (e-\y - x), ^ ~ ^ 

and 

(29) F^ oli£ (£,u) = Y, \T\W vol (Vu lT ). 

T G T(e£) : T C I? 

As announced, if ui, U2 € <5<P (£) are such that u\ = u 2 on ^T&T{eC),TcDT, then F^(C, ui) 
F £ D (C,u 2 ). 

A suitable notion to study the convergence of the functional F £ D (£) := F^(C,) is T- 
convergence. We briefly recall its definition for the unfamiliar reader. We say that a functional 
F £ : LP(D,R n ) -)■ [0,+oo] T-converges to some functional F : L p (Z),R n ) -> [0, +00] if the 
following two statements hold: 

(i) (liminf inequality) For all v G L p (D,R n ) and every sequence Vk G L p (D,R n ) which 
converges to v in L P (D, R n ) and every sequence £& — > 0, we have 

F(u) < liminf F £k (v k ); 

k— »0O 

(ii) (recovery sequence) For all v G L p (L>,R n ) and every sequence Ek — > 0, there exists a 
sequence t>fc G L P (D, R n ) which converges to t> in L P (D, R ra ) and such that 

F(u) = lim F Ek (v k ). 

The notion of T-convergence is natural for minimization problems since it ensures the conver- 
gence of minima and minimizers (recall that a set is precompact if its closure is compact): 

Lemma 4.1. [BD98, Section 7.1] Suppose that F £k : L p (D,R n ) ->■ [0,+oo] T-converges to 
some functional F : L p (D,R n ) — > [0, +00] as k — > 00. If there exists a precompact sequence 
of minimizers Uk of ' F £k in L p (D,R n ), then 

lim inf F £ , = inf F. 
fe^oo LP(D,R n ) LP(D,R n ) 

Moreover, ifuk —> ufor some u G L P (D, R n ), then u is a minimum point for F. 

For an introduction to T-convergence and its application to homogenization of integral func- 
tional we refer the reader to [Br02, BD98], whence come the following useful definitions and 
properties of integral functionals: 

Definition 4.2. We say that a function W : M nxd -> R is quasiconvex if for all A G M nxd , 

W(A) = inf J I W(A + Vu(x))dx : u G W 1,oo (D,R n )| . 

We say that it is isotropic if W(AR) = W(A) for all A G M nxa! , R G SO d . We say that it 
satisfies a standard growth condition of order 1 < p < 00 if there exists a positive constant C 
such that for all A GMT xd 

(30) C7^|A| P - C < W(A) < C(l + |A| P ). 

Let D G £>(R d ) be an open set and let W : D x M nxd -> R be a measurable function. 
We say that VF is Caratheodory if for almost every 1 6 D, A 4 W(x, A) is continuous. In 
particular, this implies that for all u G VF 1 ' 1 (D), the function x \-> VF(a?, Vu(x)) is measurable 
onD. 
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If in addition there exist 1 < p < oo and C > such that for almost every x E D, 
A i-> W(x, A) is quasiconvex and satisfies the standard growth condition (30), then the in- 
tegral functional F : W l,p (D) ->l,u4 F( u ) = Id W{x, Vu(x))dx is finite and lower- 
semicontinuous for the weak topology of W 1,P (D). 

For A E M nxd , define the function ip A : R d ->■ R n by x (->■ Ax. 

Recall that Q R := (-R,R) d for all R > 0. In [ACG], Alicandro, Cicalese and the first 
author proved the following T-convergence (or discrete homogenization) result: 

Theorem 4.2. [ACG, Theorem 5] Let < p\ < p 2 < oo. Let £ be a stationary and ergodic 
stochastic lattice in A pi n 2 . Let p > 1 and let f nn and W vo i be of class IA V and V p , respectively. 
Let D E C(R ) be an open set, and let F®(£) be the energy functional given by (27). Then the 
Junctionals F®(£) almost surely T -converge as e — > to the deterministic integral functional 

F hL : LP{D,R n ) -> [0,+oo] defined by 

nu w d , A _jf D W hom (Vu(x))dx ifueW l *(D,W>), 

I +oo otherwise, 

where Whom : M nxd — > R + is a quasiconvex function which depends only on f nn , W vo \, and 
on the stochastic lattice, and which satisfies a standard growth condition (30) of order p. Also 
for all A E M. nxd , Whom satisfies the following asymptotic homogenization formula almost 
surely: 

(32) 

W hom (A)= lim J-inf(if«(£,«) « E 6? fl (£), « = <^a on £ n Q fi \ Q fi _ 2p2 1 . 

For the link between the above result and the derivation of rubber elasticity from the statistical 
physics of interacting polymer-chains, we refer the reader to [ACG, Section 4.1], [GLTV], and 
the references therein. Note that the combination with Lemma 4.1 yields the convergence of 
minimum problems. 

Note that in (32), the lattice £ and the domain Qr correspond to £ and Dr in Theorem 2.3, 
respectively. Yet, as we shall see, the function (£, Dr) i— >• S(C, Dr) defined by 

S(C,D R ) := infJF/^Au) u E ©?«(£), u = <p A on C n D R \ D R ^ P2 X , 

is not local in the sense of Definition 2.5. 

4.2. Isotropic homogenized energy density and approximation result. Let < p\ < po < 

pi < oo, and let £ denote the random parking measure of parameter po in R d . Also, let £ R := 
£Q R be the random parking measure of parameter po in the bounded domain Q R . Let p > 1 
and let f nn and W vo \ be of class IA V and V p , respectively. 

Proposition 2.1 answers the first question of the introduction: there does exist an ergodic 
admissible stochastic lattice which is statistically isotropic, namely £: 

Theorem 4.3. Suppose n = d. Suppose there exists f nn : R + x~M. d ^ R + such that 

(33) fnn(zi,Z 2 ) = fnn(\zi\,Z 2 ) Vz 1 ,Z 2 eR d , 

and suppose also that W YO \ is isotropic. Then the energy density Whom defined by the asymptotic 
formula (32) (with £ = £ and Fe R given by (27) J is well-defined and isotropic. 

Proof. Putting £ = £ gives £ satisfying the conditions of Theorem 4.2, so the existence of 
Whom is ensured by that result. The isotropy of Wh om is a direct consequence of the isotropy 
of £ in Proposition 2. 1 and of [ACG, Theorem 9]. □ 
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We now turn to a first version of the second question: in the case of the random parking mea- 
sure, does the asymptotic formula (32) hold if £ is replaced by £, R ? This question is particularly 
relevant for numerical approximations (see [GLTV]), since the approximation of Whom on a 
bounded domain Qr also requires an approximation of the random measure £ on the domain 
Qr. The answer is positive. 

Before we state the result proper, let us provide more details on the approximation procedure. 
In Theorem 4.2, the infimum in (32) uses the Delaunay tessellation associated with the point 
process £, which is slightly nonlocal (the restriction of the Delaunay tessellation to Qr depends 
on points of £ that lie outside Qr, although in some bounded annulus around Qr due to the 
non-empty space condition). Hence, for the approximation process we need to define a local 
version of (32) through the introduction of some tessellation depending only on £ in Qr. A 
natural choice would be to consider the Delaunay tessellation associated with £ n Qr. This 
is however not suitable for the following reason: the edge lengths of this Delaunay tessellation 
are a priori only bounded by R — the Delaunay tessellation of a point set is a tessellation of 
the convex hull of the point set into d-simplices, so that there exist with positive probability 
configurations which have long edges. Arbitrarily large edge lengths are incompatible with the 
modeling of polymer chains (see [ACG, Section 4.1] and [GLTV]: in the case of unbounded 
edges, the associated functions f nn are not of class IA V , see Definition 4.1). We deal with this 
issue by considering our energy functional on a slightly smaller region than the cube Qr. 

We now introduce an approximation of Whom on Qr using £ R . The main result of this 
subsection is the following approximation result: 

Theorem 4.4. For all A G M nxd , the energy density W hom defined in (32) (with £ = 
satisfies almost surely the identity 

W hom (A) = lim -^ inf /if*- 2 " 2 (£« u) u € ©?*(£*),« = A on Q R _ 2p2 \ Q R ^ P2 ) . 
R^-oo \Qr\ u y - J 

As we shall see, Theorem 4.4 follows from Theorem 2.3 applied to (some function related 
to) the infimum in (32), seen as a function of random measures and sets. 

Remark 5. Other approximations of Whom using the random parking measure on a bounded do- 
main can be used instead of the one presented in Theorem 4.4. Let us mention two alternatives. 

First, instead of restricting the energy to the domain <3i?-2p 2 > one may consider the energy 
on the whole domain Qr provided we modify the random parking measure on Qr \ Qr-2 P2 i 
typically by taking £ R := (£ R n Qr~2p 2 ) U C^ where ( R is a suitable point set in 8Qr (say 
deterministic). The point set Q R can be chosen so that the convex envelope of £ R is Qr, and 
such that any associated Delaunay triangulation of Qr has edge lengths bounded by 4p 2 - A 
result similar to Theorem 4.4 can be proved in this case, when the boundary conditions are 
imposed on Q R only (which greatly simplifies the practical implementation, see [GLTV]). 

The second alternative is an adaptation of the popular "periodization method" in homoge- 
nization. Recall that V denotes a space-time Poisson process. For all R > we define the 
Qirperiodization of V as V^ := U keZ d (kR + Vn (Qr x R+)), and let £5 be the output of the 
graphical construction of Subsection 2.1 associated with VS in place of V. We consider a Qr- 
periodic tessellation T£ of M. d associated with £5. This allows us to define a set of Q/j-periodic 
functions on M. d : 

<5#(£#) := {u £ C(R d ,R n ) | u is ^-periodic , VT G T(^#), U| T is affine}. 
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We then let CS be a minimal connected periodic cell obtained as the union of simplices of 
T (£§.)■ Although we have not checked all the details, one should have that 



1 f C R 

W hom (A) = lim j—wi\F x *(S%,u) 

r~+oo \q r \ u y 



U G 6g($) 

such that — f5- / Vu(x)dx = A > 



almost surely. In view of the results of [GNO] on the periodization method on a simplified 
model (discrete linear elliptic equation on Z d with independent and identically distributed con- 
ductances), this approach (although it is much more delicate to implement in practice) should 
yield better convergence rates than the first two. 

The proof of Theorem 4.4 requires a series of lemmas. The first of these provides elementary 
control of energy functionals for homogeneous deformations: 

Lemma 4.5. There exists C > (depending only on d, n, pi,P2, P, fnn an d W Y0 i), such that 
for all C G A PUP2 , D G 0(R d ), all e > and all A G M nxd , F® (given by {21)) satisfies 

< F £ D ((,<p A ) < C|£>[(1 + |A[*). 

Proof. For £ G A PljP2 , no edge of T(£) has length greater than 2p2- Hence the number of 
neighbours of vertices of T(£) is bounded by some finite constant depending only on p\, p2 and 
d. This implies that the number of edges of the Delaunay tessellation T(e£) fl D is bounded 
by a constant times the volume e~ d \D\. Using Definition 4.1 and (27), this yields the desired 
result. □ 

Our next lemma shows that sequences of minimizers are precompact in L P (A, M. n ). This pre- 
compactness is a consequence of an argument a la Frechet-Kolmogorov. We give this argument 
for completeness in the following lemma since it does not explicitly appear in [ACG]. 

Lemma 4.6. Let D G C(R ) be an open set, let (v,k)k>i be a bounded sequence in L P (D, W 1 ), 
and let A G 0(D) be an open set whose closure is contained in D. If there exist an increasing 
continuous function f : M + — > R + with /(0) =0 and a sequence £k converging to zero such 
that for all h G H. with \h\ small enough, and for all k, 

(34) WhUk - u k \\LP(A,R n ) < f(\h\)+e k , 

where T^Uk : x h> Uk(x + h) (whenever it is defined), then there exists u G L p (L>,R n ) such 
that for a subsequence u& — > u in L P (A, W 1 ). 

Proof. Since the sequence u^ is bounded in L P (D, W 1 ), by the Banach-Alaoglu theorem, there 
exists u G L p (D,M. n ) such that u^ converges weakly to u for a subsequence (not relabeled). 
By the lower-semicontinuity of the norm for the weak convergence, the weak convergence of 
ThUk — Uk to ThU — u in L P (A, M. n ) for all h small enough, and by (34) 

(35) \\ThU - u\\u>(A,R") < li fc minf llr/.nfc -nfc|| L p( AR n) < f(\h\). 

Let pj be a sequence of smoothing kernels of unit mass with support in the ball B(0, 1/j) 
centred at zero and of radius 1/j. For j large enough and for all k, the convolution pj -kii^ of Uk 
with pj is in L P (A, W 1 ). In addition, by Jensen's inequality with measure pj (y)dy, and Fubini's 
theorem, 

p 

\pj * u k - u k \ p < III \u k (x-y) -u k (x)\pj(y)dy) dx 

'b(o,i/j) J 




< / \u k (x-y) -u k (x)\ p dxpj(y)dy 

lB(0,l/j) J A 
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so that by (34), 

(36) \\Pj*Uk-Uk\\LP(A,M. n ) < f(j~ 1 )+ £ k- 
Likewise, by (35) 

(37) \\pj*u-u\\lp{A,w*) < fir 1 )- 

By construction, {pj *Uk}k is an equi-continuous sequence of bounded functions in C(A). 
Indeed the functions are bounded since 

\\Pj*Uk\\L°°(A) < l!/ ? jl|L°°lkfc||Li(D,R' 1 )- 

Also for all x\, x 2 € A, 

\pj*U k (xi) -pj*U k (x 2 )\ < \X1 - X2\\\pj\\up\\u k \\Ll(D,Rn), 

so that the equi-continuity follows by the uniform boundedness of Uk in L 1 (D, MP). Hence, by 
Ascoli's theorem, for all j there exists Uj € C(A) such that pj * Uk converges up to extraction 
to Uj uniformly in A as k — > oo, and therefore in LP {A, MP). Noting that for all x G A, 



x-y)dy= j u{y)pj{x-y)dy 



lim. pj-kUk(x) = lim / u k {y)pj\^ ui<<>/ i ",,,,,,,, 

k^oo k^ooJ D J D 

by weak convergence of Uk to u, we have Uj = pj * u. Since the limit is unique, the entire 
sequence {pj * Uk}k converges to pj * u in L P (A, M n ). Hence, for all j, there exists kj such 
that for all k > kj, 

(38) \\pj*u k -pj*w|Up(a,r«) < fir 1 )- 

By the triangle inequality, for j large enough and for all k, 

\\uk-u\\LP(A,R n ) < \\uk-Pj^Uk\\LP{A^ n ) + \\Pj*Uk-pj^u\\ L p^ A ^^ + \\pj^U-u\\ LP ^ A ^ n .y 

Hence, by (36), (38), and (37), for all j large enough and all k > kj, 

\\uk -u\\lp(A,tsl"-) < 3f(j~ l ) + e k . 

This implies that u k -> u in LP [A, M n ). D 

In the proof of Theorem 4.4 we shall make use of the following T-compactness and conver- 
gence of infima results, for which we use terminology from Definition 4.2. 

Lemma 4.7. Let {CfcjfceN k> e a sequence of point sets in A Pl:P2 . Let D € 0(M d ) be an open 
set. Then for any sequence {£k}ken of positive numbers converging to zero, there exist a sub- 
sequence (not relabeled) and a Caratheodory function W : D x M. nxd — > M + such that for all 
open sets A G O(D) the Junctionals F^((k) (given by (27)) satisfy 

(39) r- lim F E A k (( k ) = F A 

k— >oo K 

where the integral functional F : L P (D, MP) — > [0, +oo] is defined by 

(40) FA{u)= U A W(x,Vu(x))dx |ftiE^(4,n 

| +oo otherwise. 

In addition, W is quasiconvex in its second variable and satisfies a standard growth condition 
(30) of order p. 

In this subsection we will only make use of this result for D = Q. 
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Proof. Lemma 4.7 is a corollary of [ACG, Theorem 3], which is the corresponding compactness 
result when the sequence (f, is a constant sequence which corresponds to the realization of a 
stochastic lattice, and for W vo \ = 0. Yet, there, the energy functionals are defined on piecewise 
constant functions (on the Voronoi tessellation associated with the point set), rather than on 
continuous and piecewise affine functions on a Delaunay tessellation. As pointed out in [ACG, 
Remark 4], the compactness result of [ACG, Theorem 3] holds as well for energy functionals 
defined on continuous and piecewise affine functions on an associated Delaunay tessellation. 
Hence, since the proof of the individual compactness result [ACG, Theorem 3] only makes 
use of deterministic arguments, it carries over to the case dealt with here provided W vo \ = 0. 
It remains to argue that one can add the volumetric term W vo \ in [ACG, Theorem 3]. Since 
this term yields a continuous contribution (it already has the form of an integral) and since 
the proof of [ACG, Theorem 3] has the same structure as for the homogenization of multiple 
integrals (application of the De Giorgi-Letta criteria, and of the integral representation results 
by Buttazzo and Dal Maso), this additional contribution can be treated in a standard way (see 
for instance [BD98, Chapter 12]). □ 

Recall notation D^ r from Section 2.2. For e > 0, r > 0, A G 0(M. d ), and ( G -A. P1>P2 , an ^ 

A € M nxd , define 

(41) © e A ' A ' r (C) := {u G 6 A (()\u = <Ph on «ni\ A 1)Cr ) U dA}. 

For all tessellations T of M, d , and all A G 0(M. d ), we define Af as the interior of the closed set 

UtgT.tcA? 1 . 

Lemma 4.8. Suppose A, A', A + G 0(M, d ) are open sets with the closure of A contained in A' 
and the closure of A' contained in A + , Let A G M nxd and let r > 2p2- Then there exists a 
constant C > such that for all ( G A pitP2 , £ G (0, 1), and u G & £ ' ,r (C)> with u extended by 
ip A on A + \ A, and all h G R with \h\ < 1/C, we have 

\\r h u - u\\ p LP{A , tRn) < C[\h\ p (l + F £ A {C,u)) + e], 

where r^ is as defined in Lemma 4.6. 

Proof. By [ACG, formula (31)] (the result can be extended to the case W vo \ > 0, since the 
associated contribution to the energy is non-negative), there is a constant C such that for all 

/iGM d with|/i| < 1/C, 

(42) \\ThU-u\\l P{A , Mn) < Cl\h\v(F £ A +((,u) + l) + e] 

for all C,e,u. By definition of the functionals and nonnegativity of f nn and W vo \, 



F £ A HCu) < F £ A ^\Cu) + F^ Ar ^(C,u) 
(43) < F e A ((,u)+F £ AAAT ^\(,u). 

Since r > 2p2, for all k, all cells of T(eC) which have a vertex in A\^ er are contained in A so 
that for u G 6 £ ' ' r (C) we have 

u= ip A on A + \A Tie0 . 

Hence by non-negativity of f nn and W vo \, 



F AAAr ^((,u) = F A+W ^(C,^a) < F/+(C,^a). 

By Lemma 4.5, the last expression is bounded uniformly in (,i £ i u - Combined with (42) and 
(43), this gives us result. □ 
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Lemma 4.9. In addition to the notation and assumptions of Lemma 4.7, let A G M. nxd and 
A G 0{D) be an open set. Let W be as in Lemma 4.7. Suppose E k is a sequence such that 
£ k > 0, linifc^oo e k = and (39) holds. Then for r > 2p2, 

lim inf {F A {( k ,u) | u G G A : A ' r {Ck)} = inf ( / W{x, A + Vu{x))dx) . 
fc ->°° ueW^ p (A) [J A J 

Proof Letr > 2p 2 . Fore > OandC G ^ Pl iP2 , define the functional F £ j4 ' A ' r (C, •) : L p (AR n ) - 
[0,+oo] andF A ' A : L p (A,R n ) ->• [0,+oo]'by 

^ £ [( " U) ~ \ +oo else; 

F A,A (u) = { F A (u) if u - <p A G W^^R"), 
\ +oo else. 
We split the proof into two steps. First we prove the T-convergence result 

(44) r- lhnF £ A k > A > r (( k ) = F A > A . 

Then we appeal to Lemma 4. 1 to prove the desired result. 

Step 1. Proof of (44). 

Let A' G 0(R d ) be an open set which contains the closure of A, and let A + G 0(R d ) be an 
open set which contains the closure of A'. 

To prove (44), we first address the r-liminf inequality. Let u k G <5 £fe ' ' r (Cfc) be a sequence 
converging to some u in L P (A, M. n ) such that 

(45) supF A > A > r (( k ,u k ) < +oo. 

A: 

Since F A (( k , •) < F £fc ' ' T (( k , •), the assumption (39) and estimate (45) imply that 

F A {u) < limm{F A {( k ,u k ) < +oo, 

fc— >oo K 

so that u G W l,p (A,M. n ). For this to imply the desired T-liminf inequality, we need to prove 
that u - if A G Wq' p (A, W n ) so that F A ' A (u) = F A (u). 

Extend u k and u by ip\ on A + \ A for all fc, so that u k — >• u in L p (^4 + ). By Lemma 4.8, and 
(45), there exists some C > such that for all fc and for all h G M. d with |/i| < 1/C, 

(46) \\ThU k -u k \\ p LP{Ai;Rn) < C(\h\ p + e k ). 

Since ThU k — u k — > r^u — u in L P (A' , W 1 ), 

|| T , ,, _ v \\p < C\h\ p 

||'/i« u 'llLP(A',R n ) — I"I 

for all /i G M d with |/t| < 1/C. From the characterization of W 1 ' p (A',W n ) by difference quo- 
tients (see for instance [Eva98, Theorem 3 Section 5.8]) we thus deduce that u G W 1,P (A' ', M n ). 
Combined with the fact that u = 92 a on A' \ A, and the continuity of the trace operator from 
W^P^A^W 1 ) to W 1 - 1 / p ' p {dA,W n ) (see for instance [Eva98, Theorem 1 Section 5.5]), this 
shows that u = <p A in W 1 ~ l / p ' p {dA, W 1 ), and therefore ii-^e W]' P (A, R n ). The r-liminf 
inequality is proved. 

To prove the existence of recovery sequences, for every u with u — ip\ G W ' p (A,M. n ) we 
have to construct a sequence u& G © £fe ' ,r (Cfc) which converges to u in L P (A, W n ) and satisfies 

F A >V) = lim F A > A > r (( k ,u k ), 

fc— >oo K 

which we may also rewrite as F A (u) = limfc^ 00 i ?A (Cfc,nfc). Using the T-convergence of 
F A (( k ), it is enough to modify a recovery sequence for u in & A k (( k ) so that it belongs to 



22 ANTOINE GLORIA AND MATHEW D. PENROSE 

©£*:' r (Cfc) an( l st iH satisfies the identity above. This modification can be achieved using De 
Giorgi's averaging method, which is a very technical argument in this discrete case. We refer 
the reader to [ACG, Proof of Proposition 3], where this issue is treated in detail. 

Step 2. Proof of the convergence of infima. 

With (44) established, we shall complete the proof by an application of Lemma 4. 1 . Let u& be 
a sequence of minimizers of F £fe ' ,r (Cfe) on ©ejj r (Cfc)- We extend up. by ip& on A + \ A for all 
k. Consider an arbitrary subsequence of the original sequence u k , with this subsequence also 
denoted u^- Since the u& are minimizers, by Lemma 4.5 we have that 

(47) supF £ A k > A ' r (( k ,u k ) < supF £ A k ' A > r (( k ,p A ) < oo. 

k k 

To check the conditions of Lemma 4.1, we shall apply Lemma 4.6. To check that we may do 
so, we need to prove that u k is bounded in L p (A + ,W l ). Since ( k G A PljP2 , and u k is affine on 
each simplex of T(efcCfc) contained in A, we claim that there exists C > such that for all k 

KE.^,, - £ £M.)r-* 

TeT(e k ( k ),TcA J 1 
- 5Z l T l SUP \ u k{x)\ P 

T&T(e k t k ),TcA * vertex of T 
< Ce d k Y, \M^kx)\ p . 

x£Ar\e k C k 

This claim follows from the following two facts: 

• the volume of any simplex in T(£k(k) is bounded by Cd{2p2£k) d (where Cd denotes 
the volume of the unit ball in dimension d), 

• the number of neighbours of any point in the tessellation is uniformly bounded (as 
proved already) so that each point x G e^Cfc H A lies in the closure of at most a fixed 
number of Delaunay cells T G T(£kCk) (independent of k). 

On the other hand, since r > 2p 2 ,fordHv G & £k ' r ((k), (v-tfA.)\A n ( ^ G W ' p (^ r ( £fcCfe )), 
and we claim that we have a discrete Poincare's inequality in 6 A k (Cfc)l~lW ' P (^T{e k ( k ))' namely 



(48) Y 4Hx)-Mx)\* 

x£Ane k ( k 



\y(y) - Vh{y) - v ( x ) + ( pajx)\ 

\y-x\ 



<c y 

{x,y} eM(ekC k ) 
\x,y) a A 

where C depends on A but not on k. The proof of (48) uses the fact, proved in [ACG, Lemma 3], 
that for any x 6 A n e^Cfc. one can find a path j(x) on the Delaunay graph from x to some 
point xq G dA-yu^) ne^Cfc with 0{e^ 1 ) steps, and moreover arrange that each Delaunay edge 
appears in at most 0{e~ k l ) of the paths j(x) for x G A n £fcCfc- Let 7(3:) = {xo, x\,...,X£,x} 
for some £ G N. Then, for any w G 6£(Ck) n W 1,p (A r(efcCh) ), 

/ M . , M V"^ , / S , M ^,V^ |W(X,-) — W(X,-_l)| 

w (x) < w(x ) +Y w ( x j )~ w ( x J-i) < c^k-^r 1 H 2 ^ 

i i x 7 J/j — 1 
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where X£ + i = x, using in addition that the edge lengths are of order e k . Since the path has 
length of order 0(e' k 1 ), Jensen's inequality yields 



w( X )\> <cJ2£k( H ? ) Z wiXi ; l)l ) 1 

j-i V \ x j x j-i\ J 



Putting w = v — tp a and summing over x G A n e k ( k , using that each Delaunay edge appears 
in at most 0{e~^ 1 ) of the paths 7(2;), yields the desired Poincare inequality (48). 

By the triangle inequality and the uniform bound \x — y\ < 2p 2 for any edge (x,y), the 
inequality (48) for v = u k yields 

e ei\ Uk{ x)T<c y, 4( luk{ tz:f )l Y+ c \M p . 



xeAne k ( k { X) y} e M{EkQ k ) 

(x,y) CA 



\y-x\ 



Using the property (25) of f nn and the definition of F £ ((k), we deduce that 

\\v,\\ p < C V r d f (r-Hv r) U k(v)-Mx) \ „ 

{x,y} £N(ek( k ) 
(x,y) CA 

< CF £ A k > A ' r (Ck,u k ) + C 

for some C > independent of k. By (47) the last expression is bounded uniformly in k. 

Using Lemma 4.8, followed by (47), we see that there are constants C, C such that for 
\h\ < 1/C and for all k we have 

\\r h u k - u k \\ p LP{A , Mn) < C[\h\ p (l + F E A h {( k ,u k )) + e k ] < C'[\h\ p + e k }. 

Hence Lemma 4.6 is applicable, showing that the sequence u k converges along some subse- 
quence to a limit in L p (>l,IR n ). Thus the closure of the original sequence u k is sequentially 
compact in L p (A,M. n ), and by the equivalence of compactness and sequential compactness in 
any metric space, this sequence is precompact. We are thus in position to apply Lemma 4.1, 
which concludes the proof. □ 

Proof of Theorem 4.4. In order to apply Theorem 2.3 to the problem under consideration, we 
need to first define the function S; we also define a family of auxiliary functions S r ,r > 0. Fix 
A e M nxd . For C 6 A PuP2 , and for e, r > 0, recall the definition (41) of & A ' K ' r {Q, and define 

U ' 

S r (Q R ,C) = inf{F?«(C,u)|uG6?«' A ' r }- 

U ' 

Let C, C' € A pitP2 with C n Q R = C n Q R . Then for any T G T(C) with T n Q R -2 P2 + 0, 
the circumsphere of T contains no point of ( in its interior by definition, so it has radius less 
than 2p2 and therefore is contained in Q R ; hence it contains no point of £' and hence T G T((') 
as well. In short the Delaunay tessellations of ( and (' coincide on Q R _2 P2 - 

Consequently, our S is local on V(Q) because changes to ( outside Q R do not affect S(Q R , Q. 
Also it can readily be deduced from Theorem 4.2 that S satisfies the averaging property on 
A Pl:P2 with respect to Q. 

It remains to prove that S is insensitive to boundary effects. Let a G (0, 1). It is enough to 
prove that for any sequence {CfcjfceN of general (pi, p 2 ) -admissible point sets and any sequence 



24 ANTOINE GLORIA AND MATHEW D. PENROSE 

of positive numbers {-R/dfcgN tending to infinity, one has 

,, m ,. \S(QR h ,Ck)-S(QR h ,jR%,Ck)\ _ 

(49) limsup — — ; = 

fc— »00 1^6-Rfcl 

where we are using notation Dr^ t from Section 2.2. In fact, we first prove this for S r , for 
r = 4p 2 - 

Set Ek = 1/Rk for all k G N. By Lemma 4.7, after extraction (we do not relabel k) there 
exists a Caratheodory function W : Q x M nxd — >■ M + such that for every open set A G 0(Q), 
(39) holds, with F A given by (40), and such that W is quasiconvex in its second variable and 
satisfies a standard growth condition (30) of order p > 1 (where p is as in Theorem 4.2). 

Next we rescale. For every open set A G 0(Q), and for all e > 0, r > 0, define 6 £ ' ,r (Cfc) 
by (41) and define the rescaling S e,r of S r by 

S^(A,Ck) = mf{F £ A (C k ,u)\ueG^ r (C k )}. 

u ' 

If R > and u G C{Qr) and if w G C(Q) is defined by v(x) = R~ 1 u(Rx), then for any 
C £ "4pi,oo we have 

F? /R (C,v)=R- d F?*(C,u), 

and moreover u G 6^ ' ,r (C) if and only if v £ &i/r V (() f° r any r > 0. Hence for any 

r > 0, 

S r (Q R ,C) = M{R d F? /R ((,v)\veeQ£ r (0} 

(50) = R d S l ^ r (Q,Ck). 

Next we claim that if (rk)keN is any sequence of numbers satisfying r& > 4/?2 and r^ = 
o(Rk) as /s — >• oo, then 



(51) limsup 



5 eh ' rfc (Q,Cfc)--RT i ^ ,a (Qfl h> Cfc) 



o. 



To prove this, first we deduce from Lemma 4.9, that for any r > 2p2 and any open set A G 
0(Q), we have 

(52) lim S £k > r {A,(k) = inf { / W(a;, A + Vu(s))da?l . 

k ~>°° ueW^ p {A) [J A J 

Let j] > 0. For all r > s > 4p 2 and all e > 0, 6?' A,r (Cfe) C 6?' A ' s (Cfc), so that by definition 
of S £ ' r and (50), for all k G N we have that 

(53) Rk d S 4p2 (Qn k Xk) = S £fc ' 4p2 (Q,C/0 < S £ ^(QXk). 

Also, by the non-negativity of the energy functions f nn and W vo \, for all k G N such that 
rk < rjRk, 

s Ek ' rk (Q,(k) < s £ ^ Rk (Q,Ck) 

(54) < F^^(( k ,cp A ) + S £ ^(Q^ v ,(k) 

provided all the simplices of T(efcCfc) which intersect Q\- v without being contained in Qi- V 
are contained in Q \ Q\-iv\ — which holds provided 2p2£fc < V since the edge lengths of the 
simplices of T(Cfc) are bounded by 2p 2 . 
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Lemma 4.5 yields the bound F®} Ql - 2,1 (( k ,Lp A ) < C\Q \ Qi-2 V \0- + |A| P ) < C'rj. Hence 
by (53) and (54), and (50) and property (52) applied to Q and to Q±- v , this implies that 



lim sup 

k— >oo 



< Cr] 



S £k ' rk (QX k )-Rk d S ip2 (Q Rk ,Ck) 

+ inf \ / W(x,A + Vu{x))dx\ - inf i / W(x, A + Vu(x))dx \ . 

We claim that the right hand side of this inequality tends to zero as r\ vanishes. It is enough 
to prove that the limsup of the difference of the infimum problems is non-positive. By the 
properties of W, u i-> f W(x, A+Vu(x))dx is continuous on W ' P {Q)- Let rj k be a sequence 

of positive numbers which tends to zero. Since C^°(Q) is dense in W ' P {Q) by definition, there 
exists a sequence u k € Cq°(Q) such that for all k € N, supp(-Ufc) C Qi-2-q k and which satisfies 

lim / W(x,A + Vu k (x))dx = inf <^ / W(x, A + Vu(x))dx > . 
k -*°°jQ ueW^ p (Q) Uq J 

Since W is non-negative and u k € W ,p (Qi~ Vk ) for all k G N, we have 

inf \ / W{x,A + Vu{x))dx, \ < / W(x, A + Vu k {x))dx, 

so that 

limsup inf < / W(x, A + Vu{x))dx > 

< inf i / W(x,A + Vu(x))dx\ , 
«eW 1 ' p (Q) iig J 

which is the claim. This proves (51). 

Now observe that Qr^rc* = Qr-r^ . The test-functions in the infimum problems defining 
S RC " +4p2 (QRXk) and o p2 (QR^Rc,C,k) coincide on Qr-r<*, since they do coincide on {T £ 
T(Cfc) | T n Qr-rc-4 P2 7^ 0}, and take the value ip A elsewhere on Qr-r<* . Hence, taking into 
account Definition 4.1, this yields for all R > 0, 

(55) S Ra+ ^(Q R ,( k ) = S 4 ^{Q R . Ra Xk) + R d - l O{R a ). 

Using (51), (50) and (55) in succession then gives us for r = 4p2 that 

?dQE k ,R?+rf n f_\ | (pd\ 

?d\ 



S r (QR k ,(k) = R a k S e ^ +r (Q,C k )+o(R a k ) 
= S R %+ r (Q Rh ,( k ) + o(R d k ) 



S r (Q Rk -R«Xk) + o(R; 



and therefore (49) holds with S replaced by S Ap2 along the subsequence, and hence along 
the original sequence too. To demonstrate (49) for S itself, note that for any R > 4/52 and 

C £ An,P2> ^ u £ ©i R ~ 2p2 ' ' (C) then extending u by tp\ to Q_r \ Qr~2p 2 gives a function 
(also denoted u) in 6>" K ' ' P2 (C), while conversely if u G qQr, a Ap2^ t h en u\q r _ 2 G 
@Qfl-2 P2 ,A 1 2p2^~) Hence there is a constant C > such that in both cases we have 

o < f?*(C,u) - Fi Qfl_2p2 (C,«) < c^- 1 , 

and hence |5(Q fl , C) - 5 4p2 (Qfi, 01 < CR 4 - 1 . Then it is clear that (49) for S follows from 

(49) for S 4p2 . 



26 ANTOINE GLORIA AND MATHEW D. PENROSE 

We are thus in position to apply Theorem 2.3, which concludes the proof of Theorem 4.4. □ 

4.3. Rubber elasticity and random parking in a bounded set. This subsection is devoted 
to a second version of the second question of the introduction, whereby we consider a more 
general domain D than the cubes considered in Section 4.2. Let D G 0(M. d ) be some fixed open 
domain. We do not focus here on the approximation of Wh om , but rather on the T-convergence 
result itself. We define the discrete model as follows. For all e > 0, we let £ 1 ' e := £ 1 / s be 
the random parking measure with parameter po > on D 1 / £ := {s~ 1 x,x G D}, let V £ be 
the associated Voronoi tessellation of M. d , and let % be the associated Delaunay tessellation of 
the convex hull of {}l £ . Although t^' £ is (po/2,2/?o)-admissible, some edges of % may be 
arbitrarily large. In order to avoid such difficulties (which are physically irrelevant anyway), we 
proceed as in the previous section and focus on Di^ £po = {x G D \ d,2(x, dD) > 4epo}- (Note 
that in the previous section D = Q so that D\^ epo = Qi : 4 £po = Qi-4 £po -) By the geometric 
argument of the proof of Theorem 4.4, £>i,4 epo has the property that every simplex T of the 
Delaunay tessellation T £ is such that if eT n Di^ £po ^ then eT C D, and therefore all the 
simplices of T £ contained in D £ -i^ po have uniformly bounded edge lengths. 

Put pi = po/2,p2 = 2po, and fix D. For each e > and ( G ^pi,p 2 (^i/e)' an( ^ °P en 

^4 ^4 

A G 0{D), we define an energy functional F £ (() := F £ ((, •) on L P (A) as: 

(56) Ff(C,u) :=F £ AnDl - 2 ^(C,u), 
where F^((, ■) is defined as in (27). We then have the following result: 

Theorem 4.10. Assume that f nn and W vo \ are of class U p and V p for some p > 1. Then 

the junctionals F £ (£ 1 ' e ) V -converge as e — > to the deterministic integral functional F^ m : 

LP(D, R n ) -)► [0, +oo] defined by 

pD M jf D W hom (Vu(x))dx ifueW^P(D,R n ), 
I +oo otherwise, 

where W^om coincides with the energy density defined by the asymptotic homogenization for- 
mula (32) of Theorem 4.2 when the random point process C is the random parking measure £ 
with parameter po in M d . 

In particular, combining Theorems 4.10 and 4.2 shows that considering the restriction of the 
random parking measure £ to Di/ £ or considering the random parking measure £ 1 ' e on Di/ £ 
yield the same thermodynamic limit for the discrete model of rubber dealt with in this section. 

To prove this result, we need the following localized version of Lemma 4.7 for point sets 
defined on bounded domains, which is proved similarly to Lemma 4.7; taking D\^ £pQ instead 
of D does not change the estimates in [ACG] since they are all set on domains compactly 
included in D. 

Lemma 4.11. Let D G 0(M. d ) be an open set, and {£k}k£N be a sequence of positive numbers 
converging to zero. Let {CfcjfceN be a sequence of (p\ , p2)-admissible point sets in F>i/ £k . Then 

there exist a subsequence (not relabeled) and a Caratheodory function W : D x M n — > M + 

^4 

such that for all open sets A G 0(D) the junctionals F £k ((k) given by (56) satisfy 

(57) r- ]imF £ A k (Ck)=F A 

where the integral functional F : L P (D, R n ) — > [0, +oo] is defined by 

(58) FA{u)= (f A W(x,Vu(x))dx ifueW^(A,R n ), 

I +oo otherwise. 
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In addition, W is quasiconvex in its second variable and satisfies a standard growth condition 
(30) of order p. 



Proof of Theorem 4.10. By Lemma 4.11, there exists almost surely an energy density W : D x 
M nxd — s- R + and a subsequence such that for all open sets A G 0{D) the functionals F £ (£ 1 / £ ) 
T-converge along the subsequence to F A : L P (A, M. n ) — > [0, +oo] defined by 



jpAf v J J A W(x,Vu(x))dx if ueW 1 ' p (A,W n ), 
I +oo otherwise. 

As in [ACG, Theorem 2] we then appeal to the characterization of non-homogeneous quasicon- 
vex functions by their minima: for almost every x G D and for all A G M. nxd , 

W(x,A) = Umj±- inf {/ W(y,A + Vv(y))\. 

P^° Np\ veWo ,p (x+Q p ,M. n ) I Jx+Q p J 

By Lemma 4.9, for all p > small enough so that x + Q 2p C D, we have 
inf { f W(y,A + Vv(y))\ 

veWg- p (x+Q p ,R n ) I Jx+Qp ) 

= Um inf |iJf f9 "(e 1/e ,«) |« € 6f +Qp (e 1/£ ),^ v?a on ^ n (x + (Q p \Q p _ 4epo ))}. 

On the one hand, as proved in Lemma 2.5, for all < a < 1 there exists an almost surely finite 
random variable Rq such that the measures £ and £ 1 / £ coincide on D^^a for all i? > i?o- Since 
for all e small enough, {e~ l y \ y G x + Q p } C D e -i E -«, there exists £o su ch that e£ and e^ 1 / 5 
coincide onx + Q p for all e < eo- Therefore, for all e < Eq, 

inf |i^ +g "(C 1/e ,t;) I v € 6r 9/, (e 1/£ ),^ = ¥>A on e^ n (x + (Q p \ Q P - 4£po ))} 

= inf |jf f9 "(e,«) | ^ G 6? +0 '(£),u = ¥>a on e£ n (x + (Q p \ Q p - 4e p ))}- 

On the other hand, as proved in Step 1 of the proof of [ACG, Theorem 2], by a suitable appli- 
cation of the subadditive ergodic theorem, 

Mm inf If? +Qp (£,v) | v G 6 £ +Qp (0,v = VA on e^n(x+(Q p \Q p ^ £po ))\ = W hom (A) 

almost surely. Hence, W(x, •) coincides with Whom for almost every x G D. In particular, by 
uniqueness of the limit, the entire sequence converges almost surely. □ 
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